Integer quantum Hall effect in a square lattice revisited 
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We investigate the phenomenon of integer quantum Hall effect in a square lattice, subjected 
to a perpendicular magnetic field, through Landauer-Biittiker formalism within the tight-binding 
framework. The oscillating nature of longitudinal resistance and near complete suppression of 
momentum relaxation processes are examined by studying the flow of charge current using Landauer- 
Keldysh prescription. Our analysis for the lattice model corroborates the finding obtained in the 
continuum model and provides a simple physical understanding. 

PACS numbers: 73.23.-b, 73.43.-f, 73.43.Cd 



I. INTRODUCTION 

The phenomenon of the integer quantum Hall effect 
(IQHE) is one of the most significant discoveries in con- 
densed matter physics^. At much low temperatures, 
typically below 1 K, this phenomenon occurs in two- 
dimensional (2D) electron systems in presence of a strong 
perpendicular magnetic field. It has been observed that 
the Hall resistance is quantized in units of h/2e 2 with a 
great accuracy, specified in parts per million, and due to 
this impressive accuracy it is utilized as a resistance stan- 
dard 2 -. Such a precise quantization of the Hall resistance 
comes from almost complete suppression of momentum 
relaxation processes in the quantum Hall regime. In a 
finite width conductor at high magnetic fields, the states 
carrying current in one direction are spatially separated 
from those carrying currents in the opposite direction, 
and as a result the overlaps between these two groups of 
states get significantly reduced which results in a suppres- 
sion of backscattering. The near complete suppression of 
backscattering between the forward and backward prop- 
agating states leads to a truly ballistic conductor even in 
the presence of impurities. The quantized nature of Hall 
resistance looks like the quantized resistance of conven- 
tional ballistic conductors, but in these conductors pre- 
cise quantization cannot be achieved since the backscat- 
tering processes are not fully removed^. 

It has also been noticed that in a Hall bar the longitudi- 
nal resistance oscillates as a function of the Fermi energy 
Ep. Naturally it may occur to our mind that when the 
Fermi energy matches with anyone of the Landau levels 
i.e., with a peak in the density of states (DOS) profile, 
the resistance should have a minimum. But the real fact 
is completely opposite to that. The resistance becomes 
minimum whenever the Fermi energy lies between two 
Landau levels where the density of states also has a min- 
imum, and, it leads to an immense curiosity and gives rise 
to the question about the current carrying states through 
the sample. It has been observed that the states those 
are situated near the edges of the sample, the so-called 
edge states, play an important role for carrying current 
when the longitudinal resistance has a minimum^. It 



is quite surprising that at the minima, the longitudinal 
resistance is almost zero and the electrons are able to 
traverse a long distance without dissipating their mo- 
mentum. It emphasizes that something special must be 
taking place which provide near complete suppression of 
the momentum relaxation processes in this regime. 

Several studies already exist in the literature which 
deal with the quantized nature of the Hall resistance, 
oscillating behavior of the longitudinal resistance, exis- 
tence of the edge states and almost complete suppression 




FIG. 1: (Color online). Four probe set-ups for determining 
(a) the Hall resistance and (b) the longitudinal resistance. 

of the backscattering processes in a Hall bar based on 
the continuum model^— . However, no rigorous effort 
has been made so far, to the best of our knowledge, to 
unravel all these features within a discrete lattice model. 
In the present work we consider a two-dimensional square 
lattice in presence of a perpendicular magnetic field and 
establish them describing the system by the tight-binding 
(TB) model. Using the Landauer-Biittiker formalism 3 - we 
find the longitudinal and transverse conductivities. On 
the other hand, to investigate the behavior of the cur- 
rent carrying states when Fermi energy lies between two 
Landau levels or coincides with a Landau level we use 
Landaucr-Keldysh prescriptio n 12 ' 13 . From this study we 
can clearly describe the momentum relaxation processes. 
The physical picture about QHE that emerges from our 
present study based on the discrete lattice model of the 
Hall bar is exactly the same as obtained in the contin- 
uum model and provides much insight to the problem 
under realistic condition. This approach could be much 
more useful in further investigation of the unconventional 
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quantum Hall effect in graphene where the Hall conduc- 
tivity is half-integer quantize d 14 ' 15 , showing plateaus at 
a xy = (Ae 2 /h)(n + 1/2), n = 0, ±1, ±2, . . . , since this 
lattice model might be adapted with more convenience 
in graphene and other topologically insulating materi- 
al a 16 ' 17 . It is our first step towards this direction. 

In what follows, we present the results. In Section 2, we 
describe the model and theoretical formulation. Section 
3 contains the numerical results and related discussions, 
and we draw our conclusions in Section 4. 



of the voltages to zero without loss of generality. Here 
we set V4 = 0. It allows us to write the currents in three 
other leads in the form of matrix equation, V=Q~ 1 I, 
where Q is a 3 x 3 matrix whose elements are expressed 
as follows: Qa — X^fc Gik and Qij = —Cry, where k runs 
from 1 to 4. The above matrix equation simplifies in the 
form: 
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II. THE MODEL AND THEORETICAL 
FORMULATION 

The quantum Hall system is shown schematically in 
Fig. [TJ where a square lattice, placed in a magnetic field 
B perpendicular to its plane, is attached to four finite 
width leads. The TB Hamiltonian for a square lattice 
with TV atomic sites in both the x and y directions reads, 



C m,n+l C m,n e m + h.C 



(1) 



Here t m ,n represents the site energy of an electron at the 
lattice site (m,n), m and n being the x and y co-ordinates 
of the site, respectively, n (c m . n ) corresponds to the 
creation (annihilation) operator of an electron at the site 
(m,n) and t is the nearest-neighbor hopping integral. The 
phase factor 9 m = 2irmBa 2 (a being the lattice spacing) 
arises due to the magnetic field B with the choice of the 
vector potential A=(0, Bx, 0). 

The side-attached leads are taken as semi-infinite 
square lattice ribbons, which are also described by TB 
Hamiltonians similar to Eq. [TJ where the hopping terms 
do not invoke the phase factor e l6m since B = in the 
leads. The lead Hamiltonians are parametrized by con- 
stant on-site potential eo and nearest-neighbor hopping 
strength to. The hopping integrals between the bound- 
ary sites of the leads and the quantum Hall system are 
parametrized by r. The number of sites in the transverse 
direction of the leads is denoted by M y . 

In order to find the Hall resistance using Landauer- 
Biittiker formalism, we connect the leads to the sample 
as prescribed in Fig.QJa) and restrict ourselves within the 
coherent transport regime 3 -. The voltage leads, namely, 
lead-2 and lead-3 are attached on the opposite sides of the 
sample. In this formalism one can treat all the leads (i.e., 
current and voltage leads) on equal footing and express 
the current in the lead-p as I p = Y^=i G pq [V p — V q ], G pq 
being the conductance coefficient related to transmission 
coefficient T pq between the lead-p and lead-q by the ex- 
pression G pq = (2e 2 / h)T pq (factor 2 accounts for the spin 
of the electrons). Here V p is the applied bias voltage in 
the lead-p. Since the currents in the various leads depend 
only on voltage differences among them, we can set one 



where, the matrix R is the inverse of Q. In our four- 
terminal set-up (Fig. HJa)) the Hall resistance is deter- 
mined from the relation, 
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On the other hand, to determine the longitudinal re- 
sistance we connect the leads to the conductor as shown 
in Fig.QJb), attaching the voltage leads lead-2 and lead-3 
on the same side of the conductor. For this set-up also 
the above expressions remain invariant, and the longitu- 
dinal resistance is given by Rl = R21 — R31 ■ Obviously, 
for this set-up the numerical values of R21 and R31 are 
different from those of the previous set-up (Fig. QJa)). 

With these Rl and Rh we can determine the longitudi- 
nal and Hall conductivities where the conductivity tensor 
a ^ is the inverse of the resistivity tensor p^ v . Mathe- 
matically it is expressed as = (pu V ) , in which p xx 
and p xy correspond to the longitudinal and Hall resistiv- 
ities, respectively. For a square lattice, p xx = p yy = Rl 
and p xy = —p yx = Rh- Therefore, we can write the final 
expressions of the conductivities asi£, 
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Next, we find the bond charge current between the 
atomic sites n and n' by the Landauer-Keldysh method 
using the relatio n 12 ' 13 , 
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Tr s [G< n ,(E)-G< n (E)}dE. (5) 
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It actually describes the flow of charges which start at 
the site n and end up at the site n' . In this expression 
G < is the conventional lesser Green's functioniSi and Vq is 
the applied potential difference between the biased leads. 
Here trace Tr s is performed in the spin Hilbcrt space, 
writing the Hamiltonian Eq. Q] including the spin of the 
electrons. 

Throughout the numerical calculations we set eo = 0, 
lattice spacing a — 1 and fix the hopping integrals (i, to 
and r) at 1. All energies arc measured in units of t and 
we choose c = e = h = 1. 
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III. NUMERICAL RESULTS AND DISCUSSION 

In the left panel of Fig. [2] we show the quantized na- 
ture of the Hall conductance together with the density of 
states for an ordered (W = 0, W measures the disorder 
strength in the sample) square lattice. The sharp peaks 
in the DOS are associated with the Landau bands those 
are produced in the presence of magnetic field B. For a 
fixed sample size, the total number of Landau bands and 
the number of states in each band strongly depend on the 
strength of the applied magnetic field. Here we choose 
the uniform magnetic field B in the form l/Q, Q being 
an integer, and it generates Q number of Landau bands 
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FIG. 2: (Color online). Left panel: Hall conductance (red 
line) as a function of Fermi energy and the electronic density 
of states (green line) for a perfect (W — 0) square lattice of 
size 20 x 20 considering B = 0.05. Right panel: Longitu- 
dinal resistance as a function of Fermi energy for the same 
parameter values. 

in the energy spectrum. This choice of the magnetic field 
makes each Landau band populated by N 2 /Q states for 
a square lattice of size NxN. Thus, for the set of param- 
eter values chosen in Fig. [3J we will get twenty Landau 
bands at the peaks of the DOS spectrum and each band 
accommodates twenty states. Since here we present the 
results only for a particular energy range, few Landau 
bands are visible. From the spectrum it is observed that 
the Hall conductance increases in discrete steps when- 
ever the Fermi energy crosses peaks in the DOS profile. 
This is what we would expect from an ordinary ballistic 
waveguide where two-terminal conductance changes in 
integer steps associated with the integer number of sub- 
bands or transverse modes at the Fermi energy Ep£. On 
the other hand, plateaus appear in the Hall conductance 
when the Fermi energy lies between two Landau bands. 
At these plateaus the Hall conductance gets the value 
(2e 2 /h)M, where M corresponds to the total number of 
Landau bands below the Fermi energy. This is exactly 
equal to the number of edge states at the Fermi energy 
and these states play the equivalent role similar to the 
transverse modes in a conventional ballistic waveguide. 

Though the quantized Hall conductance resembles to 
the quantized conductance in traditional ballistic waveg- 
uides, but the strange precise quantization can never 
be achieved in ordinary ballistic conductors since the 
backscattering processes are not fully suppressed. The 



almost complete elimination of the backscattering pro- 
cesses can only be achieved in the quantum Hall regime 
which provides surprisingly precise quantized Hall con- 
ductance. To reveal this fact, we examine the flow of 
charge current in a perfect conductor when the Fermi 
energy lies within a plateau region. The result is shown 
in the left panel of Fig. [31 where Ep is set at —2.8, and it 
predicts that the charge current flows only through the 
edges of the conductor and no current is available in the 
bulk. Very nicely we notice that the states carrying cur- 
rent from the left side of the sample to the right one are 
spatially separated from those carrying current in the 
reverse direction. As a result the overlap between the 
forward and backward propagating states is almost re- 
duced to zero leading to near complete suppression of the 
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FIG. 3: (Color online). Flow of charge current in a perfect 
(W = 0) square lattice of size 20 x 20 considering B — 0.1, 
where the left and right panels correspond to Ef = —2.8 and 
— 1.46, respectively. Left panel: Charge current flows only 
through the edges and no current is obtained in the bulk. 
Right panel: Charge current flows through the bulk. 

backscattering processes and thereby the momentum re- 
laxation process. In the continuum model these states are 
referred as the positive and negative /c-states and they are 
commonly named as the edge states. At high magnetic 
field we get practically zero overlap between these states. 
Due to this complete elimination of momentum relax- 
ation, the edge states carrying current towards the right 
side of the sample are in equilibrium with the left lead, 
while the other states carrying current in the opposite 
direction are in equilibrium with the right lead. There- 
fore, the difference in voltage measured by two voltage 
probes placed anywhere on the same side (the so-called 
longitudinal voltage) of the conductor gets zero, while a 
non-zero value of the Hall voltage is obtained when it is 
measured by two voltage probes placed anywhere on the 
opposite sides of the conductor. 

The almost entire suppression of the momentum re- 
laxation process, when Fermi energy lies in the plateau 
regions, gives rise to essentially zero longitudinal resis- 
tance exactly what we have obtained numerically (see 
the right panel of Fig. The resistance of a conductor 
is governed by the rate at which the electrons can loose 
their momentum. To relax the momentum an electron 
has to be scattered from the left side to the right side of 
the conductor through the allowed energy eigenstates in 
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the bulk of the conductor. This is practically zero (see 
the left panel of Fig. [3]) as long as the Fermi energy lies 
between two Landau levels. The non-zero value of the 
longitudinal resistance is obtained exclusively when the 
backscattering process takes place and it would be maxi- 
mum when the Fermi energy synchronizes with a Landau 
band (right panel of Fig. [5]). To articulate the authentic- 
ity, in the right panel of Fig. [31 we expose the nature of 
bond charge current in the conductor when the Fermi en- 
ergy is fixed at a Landau band {E F = —1.46). It predicts 
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FIG. 4: (Color online). Left panel: Hall conductance (red 
line) as a function of Fermi energy and the electronic density 
of states (green line) for a disordered (W — 0.5) square lattice 
of size 20 x 20 considering B = 0.05. Right panel: Longitu- 
dinal resistance as a function of Fermi energy for the same 
parameter values. 

a continuous flow of charge form one edge to the other, 
and certainly, electrons can scatter from the left to the 
right side of the conductor through its bulk. This scatter- 
ing leads to the longitudinal resistance. Thus, each time 
the Fermi energy crosses a Landau band, we get maxi- 
mum backscattering and hence a peak in the longitudinal 
resistance (see the right panel of Fig. [2]). 

We now present the results for a disordered square lat- 
tice of size 20 x 20 in Fig. 0] Disorder is introduced via a 
random distribution (width W = 0.5) of values of the on- 
site potentials (diagonal disorder), and results averaged 
over hundred disorder configurations are presented. In 
the presence of disorder sharp Landau levels are broad- 
ened (green line in the left panel of Fig. 2]), but it does 
not affect the quantized nature of the Hall conductance. 
However, in contrast to the perfect conductor, the sharp- 
ness of the Hall conductance gets reduced as observed in 
the left panel of Fig. 21 (red curve) . At the plateau regions 
the Hall conductance is again precisely quantized in units 
of 2e 2 /h and such conductance quantization, even in the 
presence of impurities, comes due to the complete spa- 
tial separation of the forward and backward propagating 
states which becomes apparent in Fig. [5] (left panel). On 
the other hand, when the Fermi energy lies on a bulk 
Landau band, electrons can easily scatter through the 
interior of the conductor, giving rise to non-zero bond 
charge currents throughout the conductor (right panel of 
Fig. [5]) . This scattering yields the longitudinal resistance 
(see the right panel of Fig. 2|)- 

It is important to note that with the increase of the 



disorder strength Landau bands get broadened more by 
virtue of disorder. As a result of disorder localized states 
are appeared in the tails of the broadened Landau bands 
and they do not contribute to the transport. But, as 
we get the non-vanishing Hall conductance, all the states 
within the Landau bands cannot be localized. At the 
band centers extended states exist which carry the Hall 
current and the amount of current which is lost due to 
the formation of localized states in the tail of a Landau 
band is exactly compensated by the remaining extended 
state d 19 ' 20 . Since the localized states within the Landau 
band do not contribute to the current the Hall plateaus 
become broadened which has been clearly explained by 
Kramer et al£^. It is also to be noted that for large 
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FIG. 5: (Color online). Flow of charge current in a disordered 
(W = 1) square lattice of size 20 x 20 considering B — 0.1, 
where the left and right panels correspond to Ef = —2.8 and 
— 1.46, respectively. Left panel: Charge current flows only 
through the edges and no current is obtained in the bulk. 
Right panel: Charge current flows through the bulk. 

enough disorder strength and for weak magnetic field, 
quantum Hall plateaus gradually vanish from the higher 
energy side. We confirm it numerically. However, it was 
not the main motivation of our present work, and the 
disappearance of the integer quantum Hall effect due to 
disorder has already been reported in the literatur o 22 ' 23 . 

Before we end this section, we would like to point out 
that the disappearance of edge states and the appear- 
ance of finite resistance associated with the backscatter- 
ing processes, as a result of movement of the Fermi level 
Ep, strongly depend on the system size, magnetic field 
and also on the impurity strength. Here we describe very 
briefly about the nature of quantum phase transition as- 
sociated with the above mentioned factors, to make the 
present communication a self contained study. The near 
zero value of the longitudinal resistance is obtained only 
when the Fermi energy lies anywhere within the plateau 
regions and in this situation the backscattering processes 
are almost suppressed which result edge currents in the 
sample. While, for all other choices of Ep, the finite value 
of longitudinal resistance is obtained i.e., the backscat- 
tering processes are dominated which provide bulk cur- 
rent. For a finite system size and fixed disorder strength, 
the appearance of distinct Landau levels, and hence the 
quantized nature of Hall plateaus, depends on the choice 
of the magnetic field. If the field is too weak, then the 
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Hall plateaus will no longer appear and the backscat- 
tcring processes are also strong enough to produce bulk 
current. On the other hand, if we fix the system size and 
choose the magnetic field to a moderate one then the ap- 
pearance of the edge currents or bulk current upon the 
movement of the Fermi level is significantly controlled 
by the impurity strength. For our system size (20 x 20) 
we examine that when W > 2, the Hall plateaus almost 
disappear and we get bulk current for any choice of Ep 
and this feature would be visible for much less value of 
W in the thermodynamic limit. Thus, we can emphasize 
that the quantum phase transition in the phenomenon of 
integer quantum Hall effect in a square lattice geometry 
strongly depends on the system size, disorder strength 
and the applied magnetic field. For a particular disorder 
strength and in presence of a moderate magnetic field, a 
functional dependence of the Fermi energy Ep with sys- 
tem size may be obtained for the crossover between the 
edge currents and bulk currents. Since at this stage it is 
really too complicated to compute it numerically we can 
try to establish this functional form in our future work, 
but, we believe, the main essence of all these features are 



quite well understood from our present study. 

IV. CONCLUSION 

In conclusion, we have investigated the phenomenon 
of integer quantum Hall effect in a discrete lattice model 
through Landaucr-Buttikcr formalism. The almost com- 
plete suppression of momentum relaxation processes 
and the oscillating nature of longitudinal resistance are 
clearly explored by examining the flow of charge current 
using Landauer-Keldysh prescription. The results pre- 
sented in this communication are worked out for absolute 
zero temperature. However, they should remain valid 
even in a certain range of finite temperatures (~ 300 K). 
This is because the broadening of the energy levels of the 
conductor due to the conductor-lead coupling is, in gen- 
eral, much larger than that of the thermal broadening^. 
Our present study about QHE based on the discrete lat- 
tice model provides much better insight to the problem 
under realistic condition. 
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